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Abstract. We study the weighted averages of resonance clusters for the 
hydrogen atom with a Stark electric field in the weak field limit. We prove 
a semiclassical Szego-type theorem for resonance clusters showing that the 
limiting distribution of the resonance shifts concentrates on the classical 
energy surface corresponding to a rescaled eigenvalue of the hydrogen atom 
Hamiltonian. This result extends Szego-type results on eigenvalue clusters 
to resonance clusters. There are two new features in this work: first, the 
study of resonance clusters requires the use of non self-adjoint operators, and 
second, the Stark perturbation is unbounded so control of the perturbation 
is achieved using localization properties of coherent states corresponding to 
hydrogen atom eigenvalues. 
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1. Introduction: Semiclassical Szego limits 



The behavior of eigenvalue clusters resulting from the perturbation of highly 
degenerate eigenvalues of elliptic operators on compact manifolds has been stud- 
ied by many researchers, notably by V. Guillemin [3] and by A. Weinstein 
The basic idea is the following. Suppose that En is an eigenvalue of an elliptic 
self-adjoint operator P with multiplicity d^, growing with N. For bounded 
perturbations Qn, satisfying ||Qjv|| — > as N — > oo, there is a cluster of nearby 
eigenvalues Enj, with the same total multiplicity as En, that tend to En as 
./V — > oo. The eigenvalue shifts unj are defined by vnj = En — Enj- The basic 
question concerns the distribution of these eigenvalue shifts as N — > oo. Since 
the eigenvalue En is increasing with N, the Hamiltonian is rescaled so that 
the eigenvalue E is independent of N. This rescaling results in a rescaled per- 
turbation and rescaled eigenvalue shifts vnj- If the rescaled perturbation Qn 
vanishes with a rate k(N), then the rescaled eigenvalue shifts un j vanish at the 
same rate. Consequently, the point measure (1/cZzv) Ylj=i ^(^ ~~ ^N,j/K(N))d\ 
should have a weak limit as N — > oo. In particular, if p € Co(R), then, roughly 
speaking, one proves that 



where A is a parameterization of the classical Hamiltonian orbits with energy 
E and \i is an invariant measure on this energy surface. The effective potential 
Q is the average of a re-scaled perturbation over one of these orbits. 

In the semiclassical context, we interpret h = 1/N as Planck's constant, 
and then this formula ([1]) is what we mean by a semiclassical Szego limit for 
the appropriately rescaled eigenvalue shifts of perturbed operator. Formula ([TJ 
expresses the weak limit of the distribution function of the eigenvalue shifts in 
terms of averages of the perturbation over corresponding classical orbits. 

The behavior of eigenvalue clusters for bounded perturbations V of the Lapla- 
cian — Agn on L 2 (S n ) were studied by Guillemin [3]. Weinstein [9] studied 
these Szego-type limits for the Laplacian on a compact manifold perturbed by 
a bounded, real-valued function V. In both cases, the semiclassical parameter 
is the index of the unperturbed eigenvalue. The integral on the right in (TJJ is 
the average of potential perturbation V over closed geodesies of the manifold. 
Brummelhuis and Uribe [2] extended these results to the study of the semiclas- 
sical Schrodinger operator — /i 2 A+V r on L 2 (M n ). The potential V > is smooth 
with Voo = liminfui^oo V(x) > 0. They studied the semiclassical behavior of 
the eigenvalue cluster near an energy < E 2 < Voo- They proved an asymp- 
totic expansion of Trp[(Hl' 2 — E)h~ l \ as h — > and related the coefficients to 
the classical flow for p 2 + V on the energy surface E. 

Uribe and Villegas-Blas [8] extended these results by considering perturba- 
tions of the hydrogen atom Hamiltonian by operators of the form e(h)Qh where 
Qh is a zero-order pseudo-differential operator uniformly bounded in h and 
e(h) = 0(h 1+s ), for 5 > 0. The main novelty comes from the fact that for a 
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fixed negative energy, there are two types of classical orbits for the Hamilton- 
ian flow on an energy surface for negative energy. There are bounded periodic 
orbits corresponding to nonzero angular momentum, and there are unbounded 
collision orbits with zero angular momentum. Uribe-Villegas [8j used Moser's 
regularization of collision orbits so that all orbits are considered periodic or- 
bits. In this regularization, all orbits correspond to geodesies on the sphere § 3 . 
Those passing through the north pole are the collision orbits. The geodesies on 
S 3 are parameterized by a certain five-dimensional set A described in Appendix 
1, section [H 

In this paper, we extend these results to resonances of the Stark hydrogen 
Hamiltonian. We prove a Szego-type result on the semiclassical behavior of 
the distribution of the resonance shifts. To explain this in more detail, let 
En{K) = —l/(2h 2 N 2 ) be an eigenvalue of the hydrogen atom Hamiltonian 
Hy(h) = -(l/2)/i 2 A - \x\- 1 i defined on L 2 (M 3 ) (see ©), with multiplicity 
d]\r = N 2 . Applying an external electric field of strength F > 0, the resulting 
Hamiltonian Hy(h,F) = Hy(h) + Fe(h)xi, called here the Stark hydrogen 
atom Hamiltonian (see ([5])), has purely absolutely continuous spectrum equal 
to the real line M. We will keep F > constant assume that e(h) vanishes as 
h — > corresponding to weak field limit. Under the perturbation by the electric 
field, the eigenvalue Ej^(h) < gives rise to a cluster of nearby resonances 
ZN,i(h, F),i = 1,...,Kn, with total algebraic multiplicity equal to dpj. We 
have $lzN,i(h, F) ~ En(K) and the imaginary part of the resonance QzN,i(h, F) 
is exponentially small in l/(hF). 

To study the semiclassical limit, we take h = 1/N as in Uribe-Villegas 
[8]. Then, the family of hydrogen atom Hamiltonians Hy(l/N,F = 0) has 
a fixed eigenvalue En (1/N) = —1/2. The Stark hydrogen atom Hamiltonian 
Hy(l/N,F) has a cluster of nearby resonances ZN t i(l/N, F) that converge to 
— 1/2 as ./V — > oo. Our main result is the following Szego-type theorem for this 
resonance cluster in the large N limit corresponding to a weak electric field. 

Theorem 1. Let F > be fixed, and let p be a function analytic in a disk 
of radius 3F about z = 0. Let e(h) = /i 6+<5 , for some 5 > 0, small, and 



take h = 1/N, with N G N. For the resonance cluster {zn,i{1/N , F)} near 
E N (1/N) = -1/2, we have 



where <pt is the Hamiltonian flow for the Kepler problem on the energy surface 
£(— 1/2) with collision orbits treated as in [8], and (pt(x(t),p(t)))i is the pro- 
jection of this flow onto the first coordinate axis x\. The measure p^ is the 
normalized Liouville measure on restricted to the energy surface S(— 1/2). 

This result parallels and extends the result of Uribe and Villegas-Blas [8j 
on eigenvalue clusters formed by bounded perturbations Qh of the hydrogen 
atom Hamiltonian. There are three new main components in this work. The 
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first is that the Stark perturbation is unbounded. The bounded perturbation 
Qh of Uribe-Villegas [8j is replaced by the unbounded perturbation Fe(h)x\ 
with small field strength as h — > 0. Control of the unbounded perturbation is 
obtained by utilizing the localization properties of coherent states of the hydro- 
gen atom Hamiltonian. The second is the fact that we work with resonances 
that appear as eigenvalues of non self-adjoint operators. Consequently, many 
estimates appearing in [8] have to be established for non self-adjoint opera- 
tors. Thirdly, we use a semiclassical result of Thomas- Villegas [7, Theorem 
4.2] to evaluate the trace of the Stark perturbation restricted to certain finite- 
dimensional subspaces (see Theorem 0) 

1.1. Contents. In section[2l we rescale the Stark hydrogen atom Hamiltonian 
using the dilation group. This establishes a countable family of rescaled Hamil- 
tonians all having a fixed eigenvalue —1/2. In section [3l we review the results of 
Herbst [I] on resonances for the Stark hydrogen atom Hamiltonian. We prove 
several important resolvent estimates necessary for our work, extending some 
estimates of Herbst [I] . The main semiclassical result is proved in section [H 
We show that the semiclassical Szego-type limit can be obtained by evaluating 
the trace of the Stark perturbation restricted to the eigenspace of the hydrogen 
atom Hamiltonian. This requires decay properties of the analytically continued 
coherent states. The final part of the proof of Theorem Q] is proved in sec- 
tion We apply a theorem of Thomas- Villegas [7] to polynomially bounded 
perturbations in order to evaluate the large N limit of the trace of the Stark 
perturbation restricted to the hydrogen atom Hamiltonian eigenspace. 
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2. Scaling 

The hydrogen atom Hamiltonian Hy{h) with the semiclassical parameter h 
acts on the Hilbert space L 2 (1R 3 ). It is the self-adjoint operator given by 



H Y (h) = ~A- X (3) 
The discrete spectrum consists of an infinite family of eigenvalues E)-(h) 

Ekih) = ^2 > k e N > W 

each eigenvalue having multiplicity k 2 . The spacing between successive eigen- 
values is 0(k~ 3 ). 

With the choice of h = 1/N and k = N, we see that E k=N (h = 1/N) = -1/2 
is in the spectra of the countable family of Hamiltonians Hy(l/N), N E N. The 
multiplicity of the eigenvalue —1/2 grows as iV 2 . 



RESONANCE CLUSTERS FOR THE HYDROGEN STARK HAMILTONIAN 



5 



The unsealed Stark hydrogen Hamiltonian is 

h 2 1 

H v (h,F) = ——A — -j— r + e(h)Fx\ 

1 \x\ 

= Hy{h)+w h (F), (5) 

where F > is the electric field strength and we have chosen the ^-direction 
for the field. We consider a parameter e(h) = h K+s , for < 5 < 1. We will 
choose K S N below. 

The dilation group D a , a > 0, a representation of the multiplicative group 
R + , has a unitary implementation of L 2 (M. d ) given by 

(D a f)(x) = a d / 2 f(ax). (6) 
We scale the Hamiltonian in ([5]) by a = h 2 : 

D h2 H v (h, F)D h - 2 =L(-1a--^ + h^Fx^j . (7) 

We call this rescaled Hamiltonian Sh(F) so that 

S h (F) = -lA-±- + h A e(h)F Xl 
2 \x\ 

= H v + h 4 e(h)Fx 1} (8) 

where we write Hy = -ffy(l) and Wh(F) = h A e{h)Fx\ is the rescaled pertur- 
bation. Note that for F = 0. the eigenvalues of S/^O) are given by E^(l) = 
-l/(2/c 2 ), with k e N. We will keep F > fixed. The effective electric field 
is h e(h)F and we can make this small be taking h — > 0, or, equivalently, with 
h = l/N, by taking N — > oo. 

3. Resonances of the Stark hydrogen Hamiltonian 

The hydrogen atom Hamiltonian Hy{h, F) with an external electric field Fx\ 
is given in ([5]). We write Ho(h,F) for the Stark Hamiltonian with V = 0. For 
i/iis section only, we take h = 1 and whie Hy(F) = Hy(l,F), for the hydrogen 
atom Stark Hamiltonian, and Hy = Hy(l,0), when the field F = 0. TTie 
Stark Hamiltonian with V = and h = 1 is denoted by Hq(F) = Hq(1,F) = 
-(1/2)A + e(/i)Fxi. and iTy = H V (1,0). For F / 0, the spectrum of Hy(F) 
is purely absolutely continuous and equal to R. We are interested in the fate 
of the negative eigenvalues of Hy when the field F is turned on. We review 
the results of Herbst [3] on the resonances associated with the F / case. 
Herbst's article deals with more general Stark Hamiltonians but we cite and 
use his results only for the hydrogen atom case of interest here for which V is 
the Coulomb potential. 

3.1. Dilation analyticity. The dilated Stark hydrogen Hamiltonian is ob- 
tained by conjugating ([5]) with the unitary dilation group D a defined in ([6]). 
We take a = e e > 0, for 9 G R. We obtain 

Hy(F,6) = D exp ^Hy(F)D exp ^_ d ) 
p—28 „—9 

= -^A-^ + e^x,. (9) 

2 \x\ 
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We are interested in extending this formula to 9 with 9$ ^ 0. There are 
two properties that need to be checked: the analyticity of the potential, and 
the analyticity of the Stark Hamiltonian. Since the potential is a Coulomb 
potential, we have V{9) = e~ 9 V, so V{9) is a type A analytic family of op- 
erators. Furthermore, as V{— A + is a compact operator, it follows that 
V{9){— A + is a compact operator- valued analytic function of 9 for 9 G C. 
In accordance with the Herbst's hypothesis [U p. 287], we may take 9q = tt/3, 
the maximum width of the strip of analyticity allowed by the purely Stark 
Hamiltonian Hq(F) = — (l/2)A + i ? xi. In order to understand the origin of the 
bound it/3, we note that for F £ M. 

H (F, 9) = -(l/2)e~ 2e A + Fe e x x = e~ 2e [-{l/2)A + Fe w Xl \. (10) 

Herbst proved results on the operator — (1/2)A + Fe 3e xi when 9 becomes com- 
plex in section II of [3] . The effective electric field Fe 3S has a nonzero imaginary 
part, necessary for Herbst's results, only if < 196*1 < tt/3. 

Herbst [4] proved the following theorem about the dilated operator Hy(F, 9). 
We consider F > fixed. 

We recall that for a closed operator with an isolated eigenvalue zq, the alge- 
braic multiplicity of the eigenvalue is defined as the dimension of the range of 
the corresponding Riesz projector. 

Theorem 2. [U Theorem III. 2] For < 96* < tt/3, the operator Hy(F,9) is 
closed on D(—A)(~)D(M Xl ). The operator family Hy(F, 9) is an analytic family 
of type A operators in 9. The spectrum on Hy(F,9) is discrete, independent of 
9, and the algebraic multiplicity of each eigenvalue is independent of 9. 

3.2. Resonances. Herbst [4 J proved that for V = 0, the closed operator Hq(F, 9), 
F / 0, has no spectrum for < \Q9\ < tt/3. As stated in Theorem [21 Herbst 
also showed, using the techniques of dilation analyticity, that for < ^s9 < tt/3, 
the non self-adjoint Hamiltonian Hy(F, 9), with V a Coulomb potential, has iso- 
lated eigenvalues with finite algebraic multiplicity. Furthermore, Herbst proved 
that these eigenvalues are connected to the eigenvalues of the F = and Q9 = 
operators. 

Theorem 3. [4, Theorem III. 3] Suppose that Eq is a negative eigenvalue of 
Hy, defined in (0) with h = 1, of multiplicity Nq. Then for F > small, 
there are exactly Nq eigenvalues, counting algebraic multiplicity, of Hy{F,9), 
as defined in (OJ) with < ^s9 < vr/3, nearby, and as F — > + , these converge 
to Eq. 

We also apply Theorems [2] and [3] to the scaled operator Sh{F) defined in 
([8]) with F of the theorem replaced by h 4 e(h)F, and take h = 1/N. For any 
fixed JV G N, we consider the resonance cluster {zNj(h, F)} of Sh(F, 9) near the 
eigenvalue — 1/(2A^ 2 ) of Hy{9). Note that the operator Sh(F) has an effective 
electric field with strength h 4 e(h)F that vanishes as h — > 0. Hence, Theorem [3] 
states that the resonances ZN,j{h,F) converge to the Af 2 -degenerate eigenvalue 
En as h — > 0. 
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3.3. Resolvent estimates. We summarize the resolvent estimates needed from 
[4]. We recall that for a closed operator A with domain D(A), the numerical 
range of A, denoted W(A), is the smallest convex set generated by {(«, Au) \ u G 
D{A)}. We let H (F) = -(1/2) A + Fx x be the Stark Hamiltonian. Following 
Herbst [3], we review the results on Stark Hamiltonians with complex electric 
fields. 

Proposition 1. [U Theorem II I] We write F = Ee^, with E, 4> G R, E ^ 0, 
anc? < |0| < 7r/3. 

(1) T/ie spectrum of Hq(F) is empty. 

(2) T/ie numerical range of Hq{F) is the half-plane 

W(H (F)) = L G C | > 3*} , (11) 

independent of E ^ 0. 

(3) 77ie resolvent is bounded 

\\(H (F) -z)- 1 1| < [^^(^(F)))]- 1 . (12) 

We now consider the dilated Stark Hamiltonian Hq(1,F,0) = Hq(F,6), as 
defined in (J9|), 

H (F,9) = -(l/2)e- 20 A + e Fx 1 . (13) 
The following operator plays an important role in the analysis: 

K(F,e,z) = v(e)(H (F,e)- z )-\ (u) 

where the dilated Coulomb potential is given by 

V{6) = fj-. (15) 
|x| 

We prove a convergence estimate for K(F,9,z) — K(Q,Q,z) with a precise 
rate of convergence as F — > (recall that h = 1 here). This estimate is possible 
since the potential is a Coulomb potential. 

We recall the basic resolvent estimates. Let Hq(9) = — (l/2)e _2e A, and 
H O (F,0) = -(l/2)e" 29 A + Fe e x 1 be the Stark Hamiltonian. For any F / 0, 
we have the following basic estimate from Proposition [TJ 

Uz-Hoi^e))' 1 ]] < l/d(z,W(H (F,6)). (16) 

Let 7at be a simple closed contour about E^ = —1/2N 2 of radius 1/(8N 3 ). For 
z G 7at, we have 

IK^-floW)" 1 !! < lM2,e' 29 M + ) = 0(iV 2 ). (17) 

The contour 7at is chosen so that it contains only one eigenvalue En of Hy(9). 
Recall that V(9)(Hq(8) + l)" 1 is a compact analytic operator valued function 
for |3?0| < rr/3. 

Proposition 2. [H Proposition III. 1] 

(1) The operator K(F,8,z) is compact and jointly analytic in (z,8) on the 
region 

{(8,z) | z G C, < \Q8\ < tt/3}. (18) 
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(2) We have the following convergence on the contour 7jv with < < 
vr/3: 

\\K(F,e,z)-K(0,9,z)\\ =0{FN i ), (19) 

as F — > 0. This convergence is uniform on the larger set {(9, z) \ d(z, W(Hq(F, 9))) > 
0, < \%9\ < vr/3}. 

We need the following lemma summarizing several key estimates on resol- 
vents. 

Lemma 1. Let z G 'Jn and < < vr/3. 

(1) \\(H (F,9)-z)-i\\ = O(N*) 

(2) \\(H (9)-z)^\\=O(N 2 ) 

(3) \\V(9)(H (9) - z)-i\\ = <D(N) 

(4) \\e- 29 Pl {H (9) - z)" 1 !! = 0(N 2 ), where Pl = -id/dx x . 

Proof. 1. The first estimate follows from the bound (112p and the fact that the 
numerical range is a half-plane located a distance 0(1/N 2 ) from the contour 

TAT- 

2. The second estimate follows similarly as the spectrum of Hq(9) is the half 
line e" 23ei IR + , see (fTBD . 

3. The third estimate requires the following bound. Let C3 denote the constant 
C 3 = (27r)- 3 / 2 (/ R3 (l + |p| 2 )- 2 d 3 py/ 2 . For any if; G # 2 (R 3 ), and for any A > 0, 
we have 

Woo<^||AV|| + C 3 A 3 / 2 |H|. (20) 

This follows from the Sobolev embedding theorem and standard estimates with 
the Fourier transform. We decompose the Coulomb potential as V = Vxb n (o) + 
V(l — Xb n (o)) = ^2 + Voo, where Xb n (o) is the characteristic function on the 
ball of radius N > centered at the origin. We have that V 2 G L 2 (R 3 ) and 

Voo e L°°(M 3 ), with ||Va|| = loI^N 1 / 2 and Halloo = 1/N. With the help of 
(f20|) . and choosing A = 1/N, we write 

||v^|| < II^IIMoc + llvyUMI 

< (^^CsllA^II + ^^ + l^HVII, (21) 

where w 3 = 4tt. Recall that (H (6) - z)- 1 : L 2 (R 3 ) -> H 2 (R 3 ) and that 
|z| = 0{N- 2 ) for z G 7at. Taking estimate with ip = (H (9) - z)" 1 ^, for 
any <j) G L 2 (R 3 ), together with estimate (2), we easily obtain estimate (3). 

4. The proof of estimate (4) follows from 

||e- 2 V(#o(0) - z)- 1 ]] < max {\ Pl (\p\ 2 - e 2e z)' 1 } = 0(N 2 ), (22) 
since z G 7at. 

□ 

We can now give the proof of Proposition [2j 
Proof. 1. As in Herbst [4], we write 

K(F, 9, z) = V(9)(-A + l)- 1 J(9, z), (23) 
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where 

J(6,z) = (-A + l)(H (F,e)-z)- 1 . (24) 

The Coulomb potential has the property that V{&){— A + l)" 1 is an analytic, 
compact operator-valued function for any 9 £ C. The quadratic estimate (|147j) 
implies that J(6>, z) is bounded. For an appropriately defined circle 7, a contour 
integral representation 



J(9,z) = (27ri)- 2 / dw d(f){w- z)'' l {<t>-9)- L J((j),w) (25) 

is used to verify that J(9, z) is analytic in the region described in the proposi- 
tion: < |90| < vr/3 and zeC. 

2. Using the resolvent formula, we write the difference on the left in (|19p as 

K(F,e,z)-K(o,e,z) 

= -Fe e x x V(e){H Q {6) - z)-\H (F, 9) - z)" 1 

+Fe 9 V(9)(H (9) - z^^O), Xl ](H (9) - z)-\H (F,9) - z)-\ 

(26) 

The commutator [Hq(9),x\] = —2e~ 2d ip\, where pi = —id/dx\. Note that 
||xiy(^)|| < e _5Re , since the potential is Coulombic. From the resolvent esti- 
mates in Lemma [H we obtain 

\\K(F,e,z)-K(6,z)\\ <Fe ne \\(H (F,e)-z)- 1 \\{A + B}, (27) 

where 

A = e- m \\{H (9) - z)- l \\ = 0(N 2 ), (28) 

and 

B = \\V(9)(H (9) - z)- 1 ]] + 2\\e~ 2d Pl (Ho(9) - z)" 1 ]] = 0(N 2 ). (29) 
Consequently, we find from (|27l) - (f29j) that 

\\K(F, 9, z) - K(0, 9, z)\\ = 0(FN 4 ). (30) 
This proves part (2) of the proposition. □ 

We recall that our F is Fh 4 e(h) so that with h = 1/N and e(h) = h K+s , part 
(2) of Proposition [2] states that uniformly for z E 7^, with |7jv| = 2n(8N 3 )~ 1 , 
we have 

\\K(F,9,z) - K(0,9,z)\\ = 0(N- K ~ S ), (31) 

as N — > 00. 



4. A SEMICLASSICAL TRACE IDENTITY FOR RESONANCE CLUSTERS 

We now return to the scaled Hamiltonian Sh(F) = Hy + Wh(F), with F > 
fixed, the perturbation W h (F) = h 4 e(h)F, with e{h) = h K+5 , and K > 6. We 
will take h = 1/N and consider iV — > 00. We need a basic trace identity relating 
the resonance shifts z N>j (l/N,F) - E N (1/N), with E N (l/N) = -1/2, to the 
eigenvalues of a reduced, finite dimensional matrix obtained from the Stark 
perturbation Wh{F). This is the main result of this section stated in Theorem 

H 
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In section [31 the dilation was written as e s . The real part of 9 does not affect 
Theorems [2] and El Consequently, we will now write the dilation as e t8 , with 
9 € M. and in the range < \9\ < ir/3. The operators are obtained by analytic 
continuation as discussed in section [3l 

As above, we fix N £ N. We study the non self-adjoint operator Sh{F,6) 
obtained from Sh{F) in flSJ) by dilation D exp ^, for 9 £ E as above: 

S h (F, 9) = D exp{ie) S h {F)D exp{ _ ie) = H V {9) + W h (F, 9), (32) 

where the dilated, scaled hydrogen atom Hamiltonian is 

e -2i6» -i0 

Hy{9) = Dexp^HyDexpt-iO) = Up ( 33 ) 

and the dilated perturbation is 

W h (F, 9) = D expm W h (F)D exp{ _ ie) = h 4 e(h)e i9 F Xl . (34) 

We write W h {F) = W h (F,0) and note that W h {F) = h 4 e(h)Fxi is self-adjoint. 

Let Il^r be the orthogonal projector for the eigenvalue = —1/(2N 2 ) 
of the scaled hydrogen atom Hamiltonian Hy defined in (jHJ). Under dilation, 
these remain eigenvalues of Hy(9). Let Pn(9) be the projector for the resonance 
cluster {zN,i(h, F)} near En of the non self-adjoint operator Sh(F, 9). We write 
for the resonance shift 

z NA {l/N, F)=E N + u N>i , u N>i G C. (35) 

The following trace estimate for the resonance shifts is the main result of this 
section. 

Theorem 4. Let vn,% be the complex resonance shifts defined in t35\) , and 
let Tjsf t i be the eigenvalues of the self-adjoint operator ITjy Wh(F)H^. For any 
m E N, and for h = 1/N , we have the following trace formula: 




for some f3 > 0. 

The proof of Theorem [3] requires two main steps. In the first, we express 
the left side of ([36]) in terms of the trace of the operator Pj^(9)(Sh(F,9) — 
EN) m P^{9). In the second step, we evaluate the trace of this operator and 
express it in terms of the finite-rank operator I\® N Wh{F)lft N . 

4.1. Step 1. A trace calculation. Since Sh(F,9) is non self-adjoint, the 
projector Pn(9) is not self-adjoint. The range of Pn(9) is a finite-dimensional 
subspace £ jy with a dimension TV 2 that is equal to the geometric multiplicity 
of the eigenvalue En of Hy. Let ZN,j(h,F), j = 1, . . . , K, with 1 < K < N 2 
be a listing of the distinct resonances that converge to En as N — > oo. Let 
Fnj(G) be the projector onto the generalized eigenspace £nj corresponding to 
the resonance ZN,j(h,F). The subspace £n admits a direct sum decomposition 
£n = ®f=\£N,ji where the finite-dimensional subspaces £ n,j, j = 1, • • • , K have 
the following properties: 
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(1) £ N . 

j — ran P^j(9) and dim£jvj — m j 

(2) Pn(9) = E?=iPnA0) 

(3) the algebraic multiplicity of the resonance znj is given by rrij 

(4) iv 2 = Ef=i^- 

(5) the projectors satisfy PN,j{Q)PN,m{6) = Sj m PN,j(9) 

(6) on the invariant subspace £nj, we have Sh(F, 9)\£jyj = ZN,jIs N +-Djvj, 
where jDjvj is nilpotent with order mj and commutes with Sh(F,9) 

(7) ker(S h (F,6)-z Njj ) m i =£ N;j . 

We refer to Kato [5, chapter III, section 6.5] for proofs of all these properties. 
The main part of the proof of step 1 is the evaluation of the trace 

Tr (p N (9)(S h (F, 9) - E N ) m P N {9)) . (37) 

Using the facts listed above, we find 

K 

Tr(p N (9)(S h (F,9)-E N ) m P N (9)) = ^Tr (p N)j (e)(S h (F, 6) - E N ) m P NJ {e) 

i=i 

K 

= Y, Tr (PnAWznA* 1 * F)-E n + D N j) m P NJ (0) 

K 

= J2(S Nd (h, F) - E N ) m TrP NJ 

K 

= Y, mj (z N>j (h,F)-E N ) m 

K 

= E m ^W (38) 
j'=i 

We also used the facts that TrP]y,j(G) = rrij, and that TrPNj{9)D]yj = 0, 
since D^,j is nilpotent and DnjPnj = Pjvj-Djvj = A/Vj- 

The distinct complex resonance shifts i/jyj are defined in (|35p . We now 
change notation and list the shifts with multiplicity included. That is, vn j is 
listed rrij times. It then follows directly from (|38H that 

( p N (9)(S h (F,9)-E N rP N (9) \ VN ,i \ m ,, q . 

4.2. Step 2. Evaluation of the trace. The second step in the proof of 
Theorem H] is to estimate the trace on the left side of (|39p . This requires that 
we replace Pn(9) by H^(9) and that we control the perturbation Wh(F, 9) 
defined in (E 



In order to replace the projector Pn(9) by 11^(9), we need some results 
relating the projector 11^(0) to Pn{@) as N — > oo, which means that h = 
l/N->0. 
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Lemma 2. For fixed 9 G M. with < 9 < tt/3, we have 

\\u° N (e)-p N (e)\\ = o(N^ K - 5 ), (40) 

for K > 3 and 5 > 0. Consequently, if P 1 - = 1 — P, we have 

(1) ||(n^)^(^(0)|| = o(iv 6 -^), 

(2) \\P^e)U%(9)\\=0(N e ~ K - s ). 

Proof. 1. We consider the contour 7jv, a circle of radius l/(8iV 3 ) > about the 
eigenvalue En- We write the projectors as contour integrals 

n%(8) - P N (8) = ^- [ dzftz-Hvie^-iz-ShfrO))- 1 ]. (41) 

Recall that S h (F,0) = Hv(9)+W~h(F, 9), and that H (h, F, 9) = -{l/2)e~ 2tf) A+ 
Wh(F,9). We define a kernel (as in section [3]) 

K h (F,9,z) = V(9)(H (h,F,9) - z)~ x . (42) 

From the second resolvent identity for Ho(h,F,9) and Sh(F,9), we obtain 

{z-ShftO))- 1 = (z-H (h,F,9)y 1 (l + K h (F,9,z))-\ (43) 

whenever the inverse on the right exists. Substituting this back into the second 
resolvent identity for Ho(h,F,9) and Sh(F,9), we obtain 

' ' ' ' ;K h (F,9,z). 



z-S h (F,9) z-H (h,F,9) z-H Q (h,F,9)l + K h (F,9,z) 

We note a similar identity for h = comparing Hy(9) = Hq{9) + V{9) with 
H {6) = -{l/2)e~ 2ie A. We let K (9,z) = V(9){H {9) - z)~ l in analogy with 
EZD for F = 0. 

' ' ' ' : K (9,z). (45) 



z - H v {9) z - H Q {9) z - H (9) 1 + K Q (9, z) 

2. We subtract (|44p from ()45|) and substitute the difference into the integral 
in flUD - Since both {z — HQ(h,F,9))~ l and (z — Hq(9))^ 1 are analytic on and 
inside 77V, their contribution to the contour integral vanishes. Consequently, 
the difference of the projections in ()4ip is equal to the contour integral 



1 

2vri 



dz 

7N 



' -K h (F,9,z) ' 



z-H (h,F,9) ay ' ' J l + K h (F,9,z) 

(46) 



' -K (9,z)- 



z-H (9) uv ' ; l + K (9,z) 

From Proposition^ part (2), we have that K^F, 9, z) and Kq(9, z) are compact 
and that Kh(F,9, z) — > Kq(0,z) in norm. We also use the fact that (z — 
Ho(h, F, #))~ 1 — > (z — Hq(9))~ 1 strongly as h — > 0. We rewrite the integrand in 
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6]) as a sum of three terms: 

1 s ,- fl ,i W | " f > |F - , -'»- wt '» lt w»,») (47) 



II 



1 1 



z-H (h,F,9) z-H o (0) 



111 s ) «■"■'> 



^•'>TTj£(M (48) 



• (49) 



i + K h (F,e,z) i + K„(e,z) 

3. We need estimates on the operators Kq(9, z) and Kh(F, 9, z) and their resol- 
vents at —1 that appear in (|4"T]). ((15]) . and (09]). We begin with the estimates 
for Kq(8,z). From part (3) of Lemma [1] we have 

\\k (9,z)\\=o(n). (so) 

From the definition of Kq(8, z), we easily find that 

(1 + K (6, z))- 1 = 1 - V(6)(H (6) + V{9) - z)-\ (51) 

Recall that the Coulomb potential V is relatively -ffo-bounded with relative 
bound less than one. So there exist constants < a < 1 and b > 0, so that for 
all u G H 2 (R 3 ), we have 

||V«|| < o||JEr u|| (52) 

Scaling by e td , with 6 G R, we find 

||V(e)u|| < o||fr (e)u|| + 6||w||. (53) 

Replacing u by R v {6)w = (H (6) + - z) -1 ^, for z G 7at and to G L 2 (R 3 ), 
we obtain 

\\V(8)R v (6)w\\ < a\\w\\ + {a\z\ + b)\\R v (9)w\\ + a\\V{9)R v (9)w\\. (54) 
Since < a < 1 and \z\ = 0(N~ 2 ), we obtain 

\\V(9)R v (e)\\<C 1 + C 2 \\R v (9)\\. (55) 
It follows from ((55]), and the fact that = 0{N 3 ) that 

IKl + i^o^^r'H^iV 3 ). (56) 

4. The second estimate we need concerns 1 + Kh(F, 9,z). We write 

l + K h (F,9,z) = [1 + K (e,z)] 

x [1 + (1 + K (9, z))-\K h (F, 9, z) - K (9, z))], (57) 

from which it follows that 

(l + KniF^z))- 1 = [l + M h (F,9,z)]- 1 [l + K (9,z)}-\ (58) 

where 

M h (F, 9, z) = (l + K (9, z))-\K h {F, 9, z) - K (9, z)). (59) 
It follows from part (2) of Proposition [2] and (j56|) that 

\\M h (F,e,z)\\ =0(N 3 ~ K - S ), (60) 
so if K > 3, this term is less than 1/2 for all N large. It follows from (I58p that 

||(1 + ir ft ( J P,0, 2 r))- 1 ]| = C(iV 3 ). (61) 



14 



P. D. HISLOP AND C. VILLEGAS-BLAS 



5. We estimate each term I, I, and III in (|47p . (|48p . and ()49|) . uniformly in 
2 G 7at, using the estimate of Proposition [2j For the first term I in (|47p . we use 
part (1) of Lemma [H part (2) of Proposition [21 and ()6ip . to obtain for i<C > 3: 

||/|| = C(iV 2 • N- R - 5 ■ N 3 ) = 0{N 5 ~ K - S ). (62) 

As for II in (|48p . we use the quadratic estimate in [H Proposition II. 4] (presented 
in appendix [8]) in order to prove the bound 

\\ Xl (H (h,F,9) - z)- 1 !! = 0(1), z G 77V . (63) 

Recalling that H (h,F,9) - H (9) = W h (F,6) = h A e(h)e id Fx u we find that 

||//|| = 0{N- 4 - K - s -N 2 -N-N 3 ) = 0{N 2 - K ~ S ). (64) 

Recalling that K > 3, we see that the term II vanishes uniformly on 77V as 
N -> 00. 

Finally, for the last term III in (149p . part (1) of Lemma [H estimate (150p . 
together with (I6ip . (156p and part (2) of Proposition O yield 

||m|| = o(n 2 ■ n ■ n 3 ■ n~ k - 5 ■ n 3 ) 

= 0(N 9 - R - 5 ). (65) 

6. The difference of the projectors on the left side of (|40p may be estimated 
from (146 p and the above estimates, recalling that |7jv| = 27r(l/(8iV 3 )): 

\\U%(9) - P N {9)\\ < h N \(\\I\\ + \\H\\ + \\IH\\) = 0(N e - K - s ), K>3. (66) 

So for K > 6, we obtain the vanishing of the difference as N — > 00. Parts (1) 
and (2) of the lemma follow from (|40p simply by writing 

(U%(9))^P N (9) = (U%(9))^(P N (9) - U%(8)), (67) 

and similarly for part (2). □ 

4.3. Step 3. Estimates on dilated coherent states. In order to con- 
trol the perturbation Wh(F,6), we need estimates on the following operators: 
U%(9)W h (F, 9)U%(9),U%(9)W h (F, 9), and the operators P N (9)(S h (F, 9)-E N )P N {9) 
and Ppf{9){Sh{F, 9) — E^)Y^ N {9). Estimates on the first two operators are given 
in Lemma El and on the second two operators are given in Lemma [H 

Heuristically, we are able to control the first two operators due to the fact that 
the matrix elements of moments of the position operator ||x|| in the eigenstates 
i> N {9) oiH v {9) satisfy 

(i, N {8),\\x\\ m 4, N {8))~N 2m . (68) 

This decay is due to the fact that the eigenstate is well localized about the 
Bohr radius and the Bohr radius scales like iV 2 . A proof of this localization 
property of the eigenfunctions is given in the Appendix 2 in section [71 This 
localization, however, is too weak to control the operator norm of the opera- 
tors IL%(8)W h (F,9)IL%(9) and U° N (8)W h (F,8). Using ([68]), we easily arrive at 
estimates of the type 

\\U N (9)W h (F,9)U° N (8)\\ = 0(N 2 e(h)h 2 ), (69) 

and the iV 2 growth will not allow control of the trace in (I36p since this is 
divided by h 2 e(h). Instead of using ([68"p . we use coherent states ^a : N that form 
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an overdetermined basis set of functions for the eigenspace of the hydrogen 
atom Hamiltonian corresponding to the eigenvalue E^. These coherent states 
were described in detail in the papers of Bander and Itzykson [1] and more 
recently in 0[8]. We recall the main points that we need here in Appendix 1, 
section [6l Recall that the dimension of the range of the projector Ti® N (9) is iV 2 
and that e{h) = h K+s , for K > 6. 

We use the following notation for operators that occur as remainder terms 
and that have bounds depending on h but uniform with respect to any other 
parameters. Let K(g(h)), for a function g(h), denote a bounded operator with 

\\K(g{h))\\=0(g{h)). (70) 

an example of a function g(h) is (e(h)h 2 ) m . We will also write K(0(h~ e )) to 
mean a bounded linear operator with [|_K"(0(/t~^))|| = 0(h ). The actual form 
of K is unimportant and may vary from line to line but a bound of the type 
rU|) or of the type 0(h~ e ) will always hold. 



Lemma 3. There exists a constant ro > 1, independent of 9 with \9\ < vr/4, so 
that for any n G N, we have 

u° N (e)D N \\\x\\ n D N2 n N (9) = u%(e)D-l\\x\\ n x M <r D N 2ii N (e) 

+U N K(O(N-°°)), (71) 

where X\\x\\<ro * s ^ e characteristic function on the set {x G R 3 | ||x|| < ro}. As 
a consequence, we have the following estimates on the perturbation restricted to 
the eigenspace of Hy{9): 

(1) \\U N (9)W h (F,9)U%(9)\\ = 0(N- K ~ 2 - S ) = 0(e(h)h 2 ) , 

(2) \\Ii%{9)W h {F,9)\\ = 0(N~ K ^- 5 ) = 0(e(h)h 2 ). 

Proof. As mentioned above, the key to controlling the perturbation is the strong 
localization property of the coherent states. Coherent states for the hydrogen 
atom are reviewed in section [6j We prove below that the dilated coherent 
states y ^ aj ]\[(e' ie x) are L 2 -valued analytic functions of 9 provided \^R9\ < 7r/4 and 
provide uniform bounds. Since, as above, Q9 pays no role in the calculations, 
we set Q9 = 0. 

1. We first prove a decay estimate for the dilated coherent states that is the 
analog of [7 ; Lemma 4.1]. We prove that there exists a constant ro > 1, inde- 
pendent of a G A, N G N, and \9\ < ir/4, so that for all n, s G N, 

lim N s [ \\x\\ n \D N 2^ aiN (e i0 x)\ 2 a?x = 0. (72) 

/ -> 00 J\\x\\>rn 



N^oo 

This estimate implies that 

% )N {e ie x)D N \ \\x\\ n x\\ x \\>r D N 2^ a)N {e ie x) d 3 x 



0(N-°°). (73) 



This estimate is uniform with respect to a, f3 G A. We turn to the proof of 
r2|) . For any b G S 2 , we will prove 

lim N s [ \x- b\ n \D N 2^ a N {e id x)\ 2 d 3 x = 0. (74) 

N ^°° J\x-b\>rn/2 
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By making appropriate choices of b we recover (|72p by a finite sum. Following 
the proof of Lemma 4.1], we first prove that for any q > 0, there exist 
constants 02,03 > 0, independent of 6, a and N, so that 



\x-b\ 2q \D N 2^ a!N {e ie x)\ 2 d 3 x 



1/2 



< 



C3Q 



\ N l/2 e c 2 N 



(75) 



(JV/240)« ' 

for any b E § 2 . We repeat the argument from Lemma 4.1] showing how 
the estimate (|75p implies (|74p . Given n and s from (|74|) . we take m E N so 
that iV — 1 < m + ra < iV and set g = (n + m)/2. Since |x • 6| > ro/2, we use 
Chebyshev's inequality and Stirling's formula for q\ to estimate (|74p from (|75p . 



AT 5 



< 



< 



< 



iE-&j>ro/2 
N s 

(ro/2)"- 

iV s 
(r /2)- 

c 5 iV s+3 



x-&| m+n |Z^2tf Q ,jv(e ie 'x)| i! tfa; 
r 3 g!iV 1 /2( A r/240)-' ? e AfC212 

jVc 4 



(ro/2)^- 1 -"" ' (76) 

where C4 = 2c2 + log 120 — 1 > 0. We choose ro > 2 so that ro/2 > e 2c4 so that 
the right side of (|76l) is bounded by 

c 6 N s+3 e- Nc \ (77) 

and this vanishes as N — > 00. This proves (I74p . 

2. To prove (j75|) . we change to the momentum variable (see (|135p ) so that 
\\(x-byD N ^ a>N (e i6 -)f 



1 ^ 



-2i0\ 



+ 1 



i0„\\JV-l 



a(N -l)(a-uj(e- w p)) 



d 3 p. 



(78) 



We next use the fact that b ■ V p generates translations in p so that for z E C, 

2 



-2i9\ 



+ 1 



a(iV- l)(a-a;(e->)) 
2 



-2i0 



(p+(z/7V)6) 2 + l 



a(iV - l)(a • w(e- ie (p + (z/N)b))"-\79) 



7V-1/ 



We need some estimates. First, in order to guarantee that the function in (|79j) 
remains in L 2 (R 3 ), we observe that if \z\/N < 1/120 and \6\ < 7r/4, there are 
finite constants < C\ , C2 so that 



2e 



2/tf 



(p + (z/AT)6) 2 + e 



< 



IRP+T 
C 2 



IIpII > 2 
IIpII < 2 - 



(80) 



This estimate is proved by estimating the absolute values of the real and imag- 
inary parts of (p + (z/N)b) 2 + e 2ie from below. So provided all other factors 
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are uniformly bounded in p, the function in (|79p is square integrable. Next, we 
prove the uniform bounds on the other factors for \z\/N < 1/2 and \6\ < tt/3. 
These conditions are less restrictive than needed for (1801) . We note that 



\a-co{e- id p)\ < VlO, (81) 

for \6\ < 7r/3 since \a\ = y/2. In order to estimate lo, we expand about p and 
write 

uj(e- ie (p + (z/N)b)) = uj{e- ie p) + V p oo(e- i9 p) ■ {z/N)e- i6 b, (82) 
for some p. It is easy to check that for \z\/N < 1/2, the gradient term satisfies 

|V p a • u(e- i6 (p + (*/i\0&))| < ci, (83) 

so that 

\a ■ oo(e- ie (p + {z/N)b))\ < Vw(l + c 2 \z\/N). (84) 
Consequently, for any N and zGCso that \z\/N < 1/2, we have 

\a-u(e- i0 (p + {z/N)b))\ N - x < c$e C2 ^, (85) 
for absolute constants cq,ci > 0. Combining these, we obtain 



2 

e- 2ie \\p\\ 2 + l 



<c 3 Ny 2 c^e c ^, 



- y "'- v r I : \ a{N -l)(a-u(e~ l9 p)) N - 1 

(86) 

since a(N - 1) ~ VN - 1. 

3. We now use Cauchy's theorem, with estimate ([86]) . in order to estimate ([75]) . 
by integrating over a path in the z-plane of radius N/240 < iV/120 about the 
origin so that estimate (I80p is valid. This gives 

\\(x-byD N2 V a>N (e ie -)\\ 

2 



< ±1 



\dz\ 

\z\=N/2m \z\ q+l 



C3 g!jVl/2 e c 2 iV 

- (A^/240)s ' 1 j 

where the finite constant C2 > is a function of cq and c\. This establishes 

(USD- 

4. It follows from estimate (|72p that for a,/3 £ A, we have 

(^.D^HxirxiWl^^*^) = 0(iV-°°), (88) 

where the error is uniform over Ax A. Of importance for us is that this estimate 
(|88p implies that the moments of the position operator in coherent states satisfy 

(v a , N (--e),D-l\\x\\ n D N 2Vp >N (--,e)) 

= (^(s0),^- 1 2 ||x|r X ||,||<ro^^(-;0)) + O(iV- 00 ). (89) 

It now follows from (|89p and the representation (|138p of the projector, suitably 
dilated, that we have the operator estimate 

u%(e)D-l\\x\\ n D N2 u%(e) = n^(0) J D iV 1 2 ||x|rxiix||<ro^n^(0)+n^(0) J R iV n^(0), 

(90) 
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where the remainder Rjy is given by 

Rn= {V a>N (-]e),D-l\\x\\ n x\\ x \\ >ro D N 2yp tN (-;6)) P a j dti(a)dfji(p), 
J A J A 

(91) 

where Pa.fi is the dyadic operator 

P a fi = \^aM e ))(^P,N^)\- ( 92 ) 

To estimate ||-Rjv||, we use estimate ([88]) and the fact that the measure \x on A 
is a probability measure. We obtain 

ll^jvll < C sup m a , N {e)\\ \\^p.N{0)\\)e- Nc \ (93) 

afieA 

where C4 > log 240 as in (|77|) . The L 2 -norms of the dilated coherent states can 
be estimates using (j80|) and (|8T1) . They satisfy the bound 

\\y a , N (0)\\ < CN 2 e N ^ 10 ^ 2 , |90| < vr/2. (94) 

Consequently, it follows from ([93]) and (fMj) that ||-Rjv|| < C~ Nc5 , for some 
c 5 > 0. Equation dHJ) then follows from flUD, and the fact that the 

measure \x on A is a probability measure 
5. We can now prove the lemma. Recall from (134)) that 

W h (F,9) = h 4 e{h)e e F Xl = h 2 e{h)D-\{F Xl )D N 2, (95) 

for N = 1/h. For part (1), we have 

\\U° N (9)W h (F,9)U° N (9)\\ < cb(fc 2 e(/i))Fro + 0(Ar°°), (96) 

for a constant r > 1. For part (2), we use \\{U° N (9)W h {F, 9))* (W h {F,9)U° N (9))\\ ■- 
\\W h (F,9)U° N (9) || 2 = \\U° N (9)W h (F,9)\\ 2 , so that estimate with n = 2 pro- 
vides the estimate. □ 

4.4. Step 4. Reduction of the perturbation. We now turn to controlling 
the operators P N (9)(S h (F,9) - E N )P N (9) and P N (9)(S h (F,9) - E N )U%(9). 
Using Lemmas [2] and [3j we can prove the analog of [8j Lemma 5] . 

Lemma 4. For any positive integer m, we have 



m,N ) 



Sh{F e [h )h ~2 EN J = (n%(9)W h (F,e)U%(9)rP N (e)+U%(9)R 

(97) 

where \\Rm,N\\ = 0(N for some f3 > 0, independent of m, and Wh(F,8) = 
e ie h 2 Fxi, with e R and < \9\ < vr/4. 

Proof. 1. To simplify notation, we suppress the 9 in the notation. We begin 
with a simple identity: 

U° N (S h - E N )P N = U° N (S h - E N )U° N P N + U° N (S h - EnXU^Pn. (98) 

We need an identity that follows from analyticity. S ince En remains an eigen- 
value of the dilated hydrogen atom Hamiltonian, Hy{9), we have 

U° N (9)(H V (9) - E N )U%(9) = 0. (99) 
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Using this f)9]), and the fact that U^HyiU ^ = 0, we have 

U° N (S h - E N )U° N = U° N W h U° N (100) 
U° N (S h -E N )(U%) ± = U° N W h (U N ) ± . (101) 
Substituting these into (J98j) , we obtain 

U%(S h - E N )P N = U%W h U° N P N + U° N W h (U%) ± P N . (102) 

2. To estimate the second term on the right in (1102j) . we use part (1) of Lemma 
[2] and part (2) of Lemma O 

\\U%W h (U%) x P N \\ < \\U%W h \\ IKK^PmW = 0(N A ~ 2K - S ). (103) 

We take the m th power of (11Q2[) and, because of (I103p . we have 

(n° N (s h - E N )p N ) m = (u N w h u° N p N ) m + n° N R m , N , (104) 

where the error term R m ^N has the form 

Rm,N = ( ^ ) llKWhKP^r-'llU^WH^^P^f, (105) 

From Lemma [2] part (1) and Lemma El we obtain 

llU^W^PMr-'llKWhiK^PNlf < C{9){e{h)h 2 ) m O{N- s ). (106) 
Consequently, R m ,N satisfies the estimate 

\\RmM\ < (e(h)h 2 ) m O(N- 5 ). (107) 

3. We next prove that 

(U° N W h U° N P N ) m = (U° N W h U° N ) m P N + U%K((e(h)h 2 ) m O(N- s )), (108) 

for all m G N by induction on m, where we use the notation K introduced 
before Lemma El We proceed by induction. Equality (|108l) is trivially true for 
m = 1. We assume it is true for m — l and verify it for m. We write 

(u° N w h n° N p N ) m = «^n^p 7V )«^<p 7V ) m - 1 

= (u%w h n% - U° N W h U N P^)[(U N W h U N ) m - 1 P N 
+U N K((h 2 e(h)) m - 1 O(N- s ))] 

= (n° N w h u° N rp N + {n%w h n%)K{{h 2 e{h)) m ~ l o{N- & )) 

^WhU^p^u^w^r^PM 

-n N W h U N P^U° N K((h 2 e(h)) m - 1 O(N- s )). (109) 

Using the estimates in Lemma [3] for LT^l^n^, we establish (|108p . 

4. We next prove a similar estimate 

(U%(S h - E N )P N ) m = U° N (S h - E N ) m P N + U° N (e(h)h 2 ) m O(N- s ). (110) 
In order to estimate the resonance term (Sh — En)Pn, we write 

(S h - E N )P N = (S h - E N )P N (P N - U%) + P N W h U%. (Ill) 
Noting that ||IIjy — Pn\\ < 1, we have from ([Hip . 

P N (S h - E N )P N = P N W h U° N (l + (U% - P N ))-\ (112) 
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From Lemma [21 we obtain the estimate 

\\(S h - E N )P N \\ <c h 2 e(h). (113) 

Given estimate (|113|) . we prove (jllUj) by induction. Assuming (jllUp for m — 1, 
we write 



(U° N (S h - E N )P N ) m - U%(S h - E N ) m P N 
= U° N (S h - E N ) [P N (U%(S h - E N )P N ) m ~ l - 

-(S h -E N ) m - l P N ] (114) 



From (I113p . we have the bound 

\\{S h -E N r~ l P N \\ < \\{{S h -E N )P N ] m - l \\ < (c h 2 e(h)) m -\ (115) 

Consequently, the norm of the left side of (|114p may be bounded above by 

||<(5 fc -^)P iV [-P i v(<) ± (^-^r' 1 ^+^n^K((/ l 2 e(/ l )) m - 1 0(iV- 5 ))] 

(H6) 

The estimate (IllOp for m now follows from this and (|113l) and (j!15|) . This 
completes the proof of Lemma H] □ 

4.5. Completion of the proof of Theorem [4]. In order to estimate the trace 
of (S h {F) - E N ) m P N on the left side in ([39]), we write 

(5 fc (F)-^rP^ = <(5 ft (F)- J B iV r^iv + (n^) ± (^(F)-^rP w . (117) 

Due to Lemma HI we have 

rrO nr /t^NttO \m /ttO \±frt / n\ B \m i_ 

N 



(S h (F)-E N rP N -(U%W h (F)U%r = (n° N ) ± (S h (F)-E N ) m P A 



+ (e(h)h 2 ) m U° N R mtN 
-(U%W h (F)U%rP^ (H8) 

We estimate the trace norm of each term on the right in (| 1 1 8 j) . 

For the first term, we use the fact that ||Pjv[|i = d^, part (1) of Lemma [2j 
and estimates on resonances in order to estimate \\(Sh(F) — E^^PnW as N^ 1 , 
and we obtain 

11(0^(^(^-^)^11! < iKn^^miP.viiaiK^CF)-^)^!! 

< d N N- a -\ (119) 

For the second term on the right in ()118|) . we have 

UemTU^^Nh < (e^^nKllillP^ivll 

< d N (e(h)h 2 ) m \\R m , N \\. (120) 

The third term is estimated as 

\\(u%w h (F)u%rpM\i < iKiixiKn^^P^rillinWll 

< d N (e(h)h 2 ) m N- a . (121) 
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Finally, we write 

Tr((S h (F)-E N ) m P N ) = Tr(U N W h (F)U N ) m 

+Tr{(S h (F) - E N ) m P N - (U N W h (F)U° N ) m } 

(122) 

with 

\Tr{(S h (F) - E N ) m P N - (U° N W h (F)U° N r}\ 

< \\(s h (F)-E N rp N -(n° N w h (F)ii%r\\ 1 

< d N {e{h)h 2 ) m N~ a . (123) 
Finally, restoring the complex parameter 8, we analyze the function 

e(0) = Tr((n%(e)W h (F, 8)U° N (8)) m ). (124) 

For 8 G R, we have, 

e(0) = rr(L» e , pW «^(F)O m J D e:Ep( _ e) ) 

= Tr((U%W h (F)U%) m ). (125) 

The function 8 — > £(0) is analytic in a neighborhood of the real axis and inde- 
pendent of 8 on the real axis, and is thus constant. Hence we can write 

Tr((S h (F, 8) - E N ) m P N (8)) = Tr((U%W h (F)U%) m ) + E m , N (8), (126) 

where E m ^{9) = 0((e(h)h 2 ) m N~ a ). This completes the proof of Theorem^) 

5. Trace estimate for the Stark perturbation of the hydrogen 

ATOM 

The next step in the proof of Theorem Q] consists of evaluating the trace on 
the right side of fl55J). Let W h (F) = h 2 Fx X - The sum on the right side of 1(55]) 
is 

note that Wh(F) = Wh(F) / (h 2 e{h)) . Since Wh(F) is a polynomially-bounded 
perturbation, we use a general result of Thomas- Villegas-Blas [7, Theorem 4.2] 
in order to evaluate the semiclassical limit of the expression on the right in 
(11271) as N -> oo. 

5.1. Polynomially-bounded perturbations. The main result of [7, Theo- 
rem 4.2] on the semiclassical limit for polynomially bounded perturbations is 
the following theorem. We slightly change notation from [7] and write h = l/N 
and k = N. 

Theorem 5. [7] Theorem 4.2] Let V be a polynomially bounded, continuous 
function on M 3 and let g : R — >• R 6e continuous. Then, we have 

c2tt 



lim -^Tr^Cn^D^-aFD^II^)) = / T V(x(t,a)) dt) d M (a). 

(128) 
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We remark that when V is a polynomial, as in our case, the proof of Theorem 
[5] is easier. In order to apply this result to (|127|) . we take g(s) = s m and 
V(x) = Fxi. Then we have from © that (D-\VD N z){x) = W h (F)(x) with 
h = l/N. 

Using the decomposition of the projector 11^ into coherent states (|138p . the 
trace on the left in (|128p . with g(s) = s m , may be expressed as an m-fold 
multiple integral with respect to ay G A of matrix elements, in the momentum 
representation, given by 

{i/N){J^ 2 K^ aiyN , (FV Pl ) J 1 / 2 ^^), (129) 

where <& a ,N is the function of § 3 defined in (I133j) , Lemma 4.3 of [7j states that 
for any 5 > 0, the matrix element in (1129j) is given by 

(i/N)(J^ 2 K^ N , J^ 2 K^ +uN ) I— F x(t, o<)i dt + OiN 6 - 1 / 2 ) 

+0(N-°°), 

(130) 

where x(t,/3)i is the first component of the vector x(t,/3) G M 3 . This vector 
and a corresponding momentum vector pit, (5), form a solution to Hamilton's 
equations for motion for the Hamiltonian h(x,p) = (l/2)p 2 — with energy 
— 1/2. The parameter (3 G A labels the Kepler orbit with energy —1/2. In the 
limit as N — > oo, the matrix elements approach zero unless cti = cti-x- This 
reduces the multiple m-fold integral to a single integral over A of the m th -power 
of the integral in (|130p . 

5.2. Conclusion of the proof of Theorem [Tl We have now proved the 
following result. For e(h) = /i 6+5 , with < S < 1 and h = l/N, the resonance 
shifts ZN,i{F, l/N) satisfy, for any m G N, the limit 

^ ( z Ntl {F,l/N)-E N {l/N) \ m f f 1 f 27r \ m 

To finish the proof of Theorem [H two steps remain to be done. First, we re- 
write the integral over A on the right in (|13ip in terms of the integral over the 
energy surface £(—1/2). Second, we show how to replace the monomial s m by 
a function p that is analytic in a fixed disk about the origin. 

As for the first task, we refer to [HJ pages 141-142]. It is noted there that the 
push-forward of the measure \i on A is the Liouville measure [il on the energy 
surface £(—1/2). Furthermore, the Kepler orbit corresponds to the Kepler flow 
4>t(x,p) on this energy surface. 

As for the second task, the function p of Theorem [T] is analytic in a disk of 
radius 3F about the origin. Since the perturbation F ■ (cj)t(x,p))i is bounded by 
2F for orbits 4>t on the energy surface £(—1/2), the estimate (|13ip guarantees 
convergence in the trace of the power series expansion for p. This concludes 
the proof of Theorem [TJ 
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6. Appendix 1: Coherent states for the hydrogen atom 

We review the construction and properties of the coherent states that form 
an overcomplete set in the eigenspace Eg of the hydrogen atom Hamiltonian 
corresponding to the eigenvalue Eg = —\/(2h?l 2 ). Let A be the five real- 
dimensional subspace of C 4 defined by 

A = {a = (pci, ai2, «3, 014) I aj G C, ||5Raj[| = ||9fa 3 -[| = 1, §Ra • Qfa = 0}. (132) 

This provides a parametrization of the co-sphere bundle S*S 3 of the three- 
sphere. There is a S'0(4)-rotationally invariant probability measure on A that 
we denote by [i. Coherent states on 8 3 have the form 

$ at e(u) = a(£)(a ■ lo) £ , weS 3 , a G A, £ G 0, 1, 2, . . . . (133) 

The coefficient a{t) ~ £ 1 / 2 is fixed by the requirement that the L 2 (S 3 )-norm of 
$> a ,N is equal to one, see [7J (2.11)]. These states are hyper-spherical harmonics. 
They are eigenstates of the spherical Laplacian — A S 3 with eigenvalue l{£ + 2). 
The entire family {^ a g(oS) \ a 6 A} is over complete and spans the eigenspace 
of — A g 3 with eigenvalue t{i + 2). We note that these states have the property 
that as iV — >• 00 they concentrate on the great circle {ui £ S 3 | \a ■ ui\ = 1} 
generated by the real and imaginary parts of a. 

In momentum space M 3 , the coherent states have the following form. The 
inverse of the stereographic projection from the three sphere S 3 to M 3 is the 
mapping p£l 3 -> co(p) G § 3 defined by 

For any a £ A, we define 

KM = <*■ - !)^ /2 ( £ 2| H | 2 + 1 )'( Q • u{lp)Y-\ P G M 3 . (135) 

These functions are in L 2 (M 3 ). Their Fourier transforms are eigenfunctions of 
Hy with eigenvalue Eg [U section II. B]. They form an overdetermined basis of 
the £ 2 -dimensional eigenspace Eg C L 2 (M 3 ) in the momentum space representa- 
tion. 

For 9 G R, we scale these momentum space functions by e e to obtain 
^(e e p) = a(£ - IK 3 / 2 (« • u(£e 9 p)y-\ p G R 3 . (136) 

These functions are analytic L 2 -valued functions for \$s9\ < tt/2. 



2p 



\\P 



| 2 + r 
| 2 -i 
| 2 + i' 



3 = 1,2,3 



(134) 
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The configuration space coherent states are obtained by the inverse Fourier 
transform: 

- (^) V^)>- 

(137) 

These form an overdetermined basis of normalized (but not orthogonal) L 2 - 
functions for &i C L 2 (R 3 ) in the configuration space picture. Let fi be the 
probability measure on A. The orthogonal projector LT^ onto the eigenspace £i 
may be written as 

U° e = £ 2 [ |*«^)(* a ,/| dfji(a). (138) 



>A 

For # G M, the dilated coherent states are 
^a,i{x-0) = D exm ^ a/ (x) = e 3d/2 ^ a/ (e d x) 

-36/2 

d*p e™* * a 4e-"p) 



(2tt) 3 / 2 



From the comments after (|136p . it follows that ^ a ^{x\9) has an analytic con- 
tinuation as an L 2 -valued function of 9 for < tt/2. We also consider the 
dilated projector H®(9) = D exp rg\Il^D exp fgy We extend these dilated opera- 
tors to 9 £ C provided |9f0| < ir/2. As mentioned above, the coherent states 
concentrate on the Kepler orbit as £ — > oo. Lemma [3] is a consequence of this 
fact. 

7. Appendix 2: Moments of the hydrogen atom eigenfunctions 

We give the brief proof of the localization property of the eigenstates of the 
hydrogen atom Hamiltonian discussed in section I4.3I The eigenfunctions of 
H v may be written as VWmO*, 0, 4>) = e~ r/n F n i(r)Yi m (9, <p), with n= 1,2,..., 
the principal quantum number. The functions Y^ m are the spherical har- 
monics {£, m) labeling the angular momentum so that < £ < n — 1 and 
—£< m < £. The radial component F n i(r) = A n i(2r/n)^L 2 ^'}_ 1 (2r/n), where 
^-£-i( r ) is the associated Laguerre polynomial of degree n — £ — 1. The 
normalization constant is A n £ = (2/n 2 )^/ (n — t — l)!/[(n + £)l] 3 . The dilated 
eigenfunctions ip n e m (9)(r) = (D expW ^ nem )(r) = e 39 / 2 ip nem (e e r) are analytic 
L 2 -valued functions in the strip \^s9\ < tt/2. Consequently, the function f(9) = 
(VWm(#), \\x\\ k tp n £' m '(9)) is analytic on the strip \Ss9\ < tt/2. For 9 G R, we have 
f{9) = e~ ke f(0). The function f(9) = e~ ke f(0) is an entire analytic function. 
By the identity principle for analytic functions, we have f{9) = f(9) on the 
strip \Q9\ < tt/2. 
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We first want to estimate the large n behavior of the expectation of powers 
of the position operator in the dilated eigenstates VWm as in (|68j) . Since f(9) = 
e~ me /(0), this reduces to an estimate on moments of the position operator. For 
any k > —21 — 1, these moments are well-known and computed in section 3, 
Appendix B of Messiah [6] . Let 

{\\x\\ k )n£m = (VWmW, INI* VW'm' (#)) ■ (140) 

Note that because of the spherical symmetry, the quantum numbers (£, m) must 
be the same in (|140p . Properties of the Laguerre polynomials leads to 

<IMIW = ^[3n 2 -£(£ + l)], (141) 

and a recursion formula for higher moments k > 2: 
9k -4- 1 n^k 

{\\x\\ k ) nem = n 2 f^(IM| fc - 1 >^ m - i ^ I y[(2^+l) 2 -A ; 2 ](||xf- 2 ) n , m . (142) 

It is easy to verify from (|14ip and (|142j) that 

(Vw m JxfVw(#)> = 0(n 2k ), (143) 

since < I < n — 1 (see [H section B.3]), verifying ([68|) . 

Secondly, we want to compute the matrix elements of the Stark perturbation. 
It is easier to do the calculation with the electric field in the three direction 
#3 because of the properties of the usual spherical harmonics. The result is 
independent of the field direction so we assume this here. Recalling the identity 
principle for analytic functions used above, the matrix element of the Stark 
perturbation is 

(^nem{0),W h (F,9)ip ne , m ,(8)) = e(h)h 2 e e F(i; nem (9),x 3 ^ M (9)) 

= e(/i)/i 2 F(T/w m ,x 3 Vw'm'>, |9?0|<7r/2. 

(144) 

We write X3 = rcosQ, where 9 is the angle with the X3-axis. By the Cauchy- 
Schwarz inequality and the first moment estimate (|14ip . we have 

\(ipnem,X3i/j n e'm'}\ < \\r 1/2 ipni'm'\\ \\r 1/2 ipnim\\ = 0(n 2 ), (145) 
since < / < n — 1 (see (6J section B.3].) 

8. Appendix 3: The quadratic estimate 

We restate the quadratic estimate appearing in Proposition II. 4 of [4]. In the 
notation of that paper, let h(a) = —A + ax\ where \Qa\ > 0. Let 9 = arga 
and define two constants 

c(a) = (3/2)(l - |cos0|) 3/2 |sin#| 3/2 H 4/3 , /3(a) = (1 - |cos0|)/2. (146) 
then, for all tp E S(R 3 ), we have 

\\h(a)n 2 + c(a)\\n 2 > /?(«)(||A^|| 2 + ||x^|| 2 ). (147) 
In our case, we have 

H (h, F, 9) = e- 2lt) [-A + e 3ie /i 4 e(/i)F 2 ; 1 ] = e~ 2ie /i(a), (148) 
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where a = e 3j6l /i 4 e(/i)-F. From ()146|) . we find that 

C(fi) = 0(iV-( 4 /3)(4+K+5) ); /?( ~ ) = Q{l y 



(149) 



Consequently, we have the following bound for z £ 7 at: 



x 1 (H (h,F,9)-z) 



i 



0(1). 



(150) 
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